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ABSTRACT It is known that lipid membranes become permeable in their melting regime. In microscopic con- 
ductance measurements on black hpid membranes one finds that conduction takes place via quantized events closely 
resembling those reported for protein ion channels. Here, we present data of ion currents through black lipid membranes 
in the presence and absence of the anesthetics octanol and ethanol, and compare them to a statistical thermodynamics 
model using parameters that are obtained from experimental calorimetric data. The conductance steps in pure lipid 
membrane suggest aqueous pores with the size of approximately one lipid cross-section. We model the permeabiUty by 
assuming empty sites of the size of one lipid. We find that pore formation in the melting transition regime is facilitated 
by the increase of the lateral compressibility that expresses itself in the area fluctuations. Thus, pore formation is related 
to critical opalescence in two dimensions. Anesthetics alter the permeability by affecting the thermodynamic state of the 
membrane and by shifting the heat capacity profiles. 
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1 Introduction 

In biology, lipid membranes are typically considered as insu- 
lators, which is an important assumption for models related to 
nerve pulse propagation (e.g., the Hodgkin-Huxley model 1 1 1) 
and the conduction of ions through proteins. It is known, how- 
ever, that lipid membranes become permeable in the proximity 
of membrane melting transitions ||2l|3]|4]|5], which creates sig- 
nificant conceptual problems for such models because biolog- 
ical membranes are in fact often close to such transitions un- 
der physiological conditions (|6l. The increased permeability of 
membranes in transitions has been attributed to the large fluc- 
tuations in membrane area generally related to transitions, and 
to the appearance of domain boundaries |2 , 7|. The area fluctu- 
ations are proportional to the lateral compressibility 1 8 1 and the 
compressibility is therefore high in the transition regime. Thus, 
the work necessary to create conducting membrane defects un- 
der these conditions is small |3 1. 

Interestingly, it has been found in microscopic measurements 
using black lipid membranes or using patch-clamp that the con- 
ductance of ions through pure lipid membranes occurs in quan- 
tized steps, suggesting pores of well-defined size ll9l [T0lfm[T2l 
|T3][l4ll21. These events display a similar conductance and typi- 
cal opening and closing time scales as the conduction events re- 
ported for protein ion channels [ 15!|. In a recent paper we found 
conductance steps suggesting aqueous pores of 0.7 nm, very 



similar to the size proposed for protein ion channels f5\. Pores 
of this size have also been proposed on the basis of molecular 
dynamics simulations [16|. The obvious question arises how 
the lipid conduction events can be distinguished from those of 
proteins. In particular, while lipid membranes can be inves- 
tigated in the absence of proteins, channel proteins cannot be 
investigated in the absence of lipid membranes into which they 
are embedded. A patch pipette has a diameter of about 1 fim 
while the diameter of a typical membrane protein is of the or- 
der of 5 nm. The patch tip area is therefore about 40000 times 
larger than the protein, which renders it technically impossible 
to measure currents through proteins without simultaneously 
measuring the events in the lipid membrane, too. 
Recently, using fluorescence correlation spectroscopy we 
showed empirically that the permeability of membranes is with- 
in experimentally accuracy proportional to the heat capacity in 
the transition range [5|. Thus, everything that shifts the heat 
capacity profile of membrane melting will also alter the mem- 
brane permeability. This includes changes in temperature, pres- 
sure and pH, but also the presence of anesthetics. The latter 
molecules have been shown to lower membrane transitions ac- 
cording a freezing point depression law, completely indepen- 
dent of the nature of the anesthetics drug [17]. When mea- 
suring in the transition regime, anesthetics therefore lower the 
permeability, and quantized current events are blocked [5] 
In this paper we show experimental data on quantized currents 
through black lipid membranes (BLJVls) and the effect of anes- 
thetics on the conduction. We rationalize our findings with a 
statistical thermodynamics model that takes the area fluctua- 
tions of the membranes into account. Statistical thermodynam- 
ics models (e.g., 2 state Ising models or 10 state Pink models) 
in the past have been successfully used to describe cooperative 
melting events in membranes DSl|T9l|20l|2l]|22l|23]|2l. In 
transitions, the cooperative fluctuations lead to domain forma- 
tion on scales much larger than the individual lipid. On these 
scales, the lack of molecular detail is not important. For this 
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reason,. Here, we allow for the formation of pores facilitated 
by the fluctuations in area. We demonstrate that this models 
results in channel statistics similar to the BLM experiment and 
to explain the effect of anesthetics. 

2 Methods 
2.1 Materials 

Chemicals. Decane and 1-octanol were purchased from Fluka 
Chemie AG (Deisenhofen, Germany), N-hexadecane, chloro- 
form and ethanol were obtained from Merck (Hohenbrunn, Ger- 
many), n-pentane was provided by BDH (Poole, UK) and potas- 
sium chloride from J.T. Baker Analyzed (Deventer, Holland). 
l,2-dipalmitoyl-sn-glycero-3- phosphocholine (DPPC) and 1,2- 
dioleoyl-sn-glycero-3-phosphocholine (DOPC) were purchased 
from Avanti Polar Lipids (Birmingham, AL) and used without 
further purification. For all experiments, MilliQ water (18.1 
Mil) was used. 

Calorimetry. DSC experiments were performed on a VP-DSC 
Calorimeter (MicroCal, Nort-hampton, MA, USA) with a scan 
rate of 5°C/h. The lipid samples were prepared by pre-dissol- 
ving in chloroform, which was first dried under a nitrogen 
stream and the kept under high vacuum overnight. The dried 
lipid mixtures were dispersed in MilliQ water to a final con- 
centration of 20mM. The buffer used was the same as for the 
BLM experiments. Before filling the calorimeter, the solutions 
were degassed for 10 minutes. 

Black lipid membranes. Planar bilayers were formed over 
a round aperture of «80 /im radius in a Teflon film of 25 /im 
thickness dividing two compartments of a teflon chamber em- 
bedded in a brass block that could be heated by a circulating 
water bath. The aperture in the Teflon film was prepainted 
with 5% hexadecane in pentane. The BLMs were painted with 
DOPC:DPPC 2:1 lipid solutions in decane/chloroform/metha- 
nol 7:2:1 and formed following the method described by Mon- 
tal and Miiller ll25l . The two compartments of the Teflon block 
were filled with unbuffered 150mM KCl (pH 6.5). Lipid solu- 
tion (25mg/ml) was spread on the buffer surface in each com- 
partment (approx. 3 /il on each side). Ag/AgCl electrodes were 
placed into both compartments of the chamber After 15-30 
minutes to allow for the evaporation of the solvent, the wa- 
ter level of the compartments was lowered and raised several 
times until a bilayer was formed over the hole. The formation 
of a BLMs was controlled visually and by capacitance mea- 
surements (with triangular lOOmV voltage input pulse). For 
the experiments with octanol, 15% v/v octanol in methanol so- 
lution was prepared and added symmetrically on both sides of 
the Teflon chamber. To calculate the concentration of octanol 
in the membrane we assumed an octanol partition coefficient 
between lipid membrane and water of 200 |26|. In the exper- 
iments with ethanol we used 60 /il and 120/il of 99% ethanol 
for a buffer volume of the two compartments of 2ml. Assuming 
a ethanol partition coefficient of 0.48 between membranes and 



water fTl] this corresponds to 9.5 mol% and 19 mol % ethanol 
in the membrane. 

Conductance measurements were performed on an Axopatch 
200B amplifier in voltage clamp mode connected to a DigiData 
1200 digitizer (both Molecular Devices, Sunnyvale, CA, USA). 
Current traces were filtered with IkHz low-pass Bessel filter 
and recorded with Clampex 9.2 software (Axon Instruments) 
on the hard drive of the computer using an AD converter with 
a time resolution of 0.1ms. The data was further analyzed with 
Clampfit 9.2 and low-pass filtered with Bessel (8-pole) filter at 
a cut-off frequency of 300Hz. Temperature was controlled by a 
HAAKE DC30 K20 (Waltham, USA) waterbath and a thermo- 
couple (WSE, Thermocoax). 



2.2 Simulating pore formation 

Here we introduce a statistical mechanical model for the de- 
scription of pore formation close to lipid melting transitions in 
a system with one species of lipids and a second species of 
molecules (anesthetics), based on previous work by ll24ll . In 
this work each lipid could be in either gel or fluid state. This 
model formally corresponds to an Ising model in a field. In the 
present work we additionally allowed for the formation of pores 
in the membrane by moving lipids from randomly chosen lat- 
tice sites to the edge of the computer matrix. For a given tem- 
perature we evaluate the mean number of gel and fluid lipids 
and the fluctuations in their number, and we determine the num- 
ber of spontaneously formed pores. 

A numerical evaluation by means of Monte Carlo simulations 
was used. To do this one first has to define the various Monte 
Carlo steps. Whether a given step is accepted or not, depends 
on difference in the free energy of the old and the new sys- 
tem configuration as given by eq. (1). At its most basic (i.e. 
a one-component system without anesthetics and no pore for- 
mation) only one type of step is needed, namely one that al- 
lows the lipids to change their state (a "melting step"). As 
we are using the two-state Ising model as a basis, each lipid 
is assumed to be in either an ordered (gel) or disordered (fluid) 
state t2n l24l . For the more advanced models (multiple lipid 
species, anesthetics, pore formation, etc.) it is also necessary 
to include the possibility for any two particles to swap posi- 
tions (i.e., a "diffusion step", fSS', '29]), as well as to include 
a Monte Carlo step which allows for the creation/sealing of a 
pore. The basic idea behind the pore forming step is the ob- 
servation that a lipid changes area by approximately 25% when 
it goes from the fluid to the gel state. So if three neighboring 
fluid lipids simultaneously change states, one can create a pore 
which has an area that is equal to a gel state lipid. Thus this 
step both conserves the number of lipids and also the area as 

Abe fore = SAfiuid^ Aafter = SAgel + Apore- In this kind of 

system a pore has automatically the cross-section of one lipid, 
which is what is observed in experiments [5J and molecular dy- 
namics simulations lfT6l . 
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Figure 1: Left and center: The pore formation step. Green represents fluid state lipids, while the gels are red. The yellow 
particle can be either a lipid or an anesthetic molecule. In this Monte Carlo step three lipids simultaneously change state while 
the yellow particle is moved to or from the edge of the lattice. This means that the pore formation step increases the lattice size 
by one. During this step the overall area of the matrix stays constant (see text). Right: The dark and light grey regions (bottom 
and top half) indicate the boundaries of repeating lattices. The definition of the nearest neighbors of the particle moved from 
the site of the pore to the edge. The red lines indicate the definition of nearest neighbors at the edge. Due to the odd shape of 
the lattice the blue lattice sites display a different definition of nearest neighborhood (indicated by the dotted red line). Still all 
lattice sites have 6 nearest neighbors. 



First pore formation is described. This step can be applied if 
a matrix site (that is not already a pore) is randomly picked: 
first the surroundings (nearest neighbors) is checked for lipids 
in the fluid state. If there are three or more, a pore is temporar- 
ily created by changing the three fluids into gels, and moving 
the "central particle" to the "end" of the lattice, increasing the 
number of lattice sites by one (see Fig. [T] left and center). In 
order to obey detailed balance, there must also be a possibil- 
ity to seal a pore. This can happen when a pore is picked, and 
it has three or more neighbors which are gel state lipids. If 
so, the pore is temporarily sealed by changing the three gels 
to fluids, and moving the "end particle" into the pore, thus re- 
ducing the number of lattice sites by one (see Fig.[T]i. Finally, 
the proposed change is accepted or rejected using the Glauber 
algorithm BOJ . We used periodic boundary conditions using a 
special definition of the nearest neighbor lipid environment on 
the edge (see Fig.[T| right). 
Model parameters 

There are three parameters which need to be determined if one 
wants to describe the melting transition of the most basic one- 
component two-state lipid system, namely the melting enthalpy, 
A_ff , the melting entropy, AS", and the cooperativity parameter, 
uofg. The latter parameter describes the interaction between a 
gel and a fluid lipid and determines the half width of the melt- 
ing transition. For the model additionally containing anesthet- 
ics and pores, another five interaction parameters are needed, 
namely Uaf (anesthetic/fluid), ojag (anesthetic/gel), LOap (anes- 
thetic/pore), ujpg (pore/gel), and Wpf (pore/fluid) ( ll23l ). 
With these parameter given the changes in free energy can be 
calculated as 

AG= AH ~TAS + ANfgLUfg 

+ ANafCUaf + ANagLOag (1) 
+ ANpfUJpf + ANpgUJpg + ANapUJap 

Fortunately, the first three parameters can be obtained from 



calorimetric experiments ( ll24l ). while the remaining three will 
need to be guesstimated from permeation measurements. Specif- 
ically, one can obtain the melting enthalpy and entropy of the 
lipids from the heat capacity profile of the one-component lipid 
system: the enthalpy change of the transition is obtained by in- 
tegrating over the excess heat capacity in the relevant tempera- 
ture interval. The entropy change can then easily be determined 
via AS" = AH/Tm, as the transition midpoint, T„i, can be read 
off the heat capacity profile directly. The cooperativity parame- 
ter, ojfg, relates directly to the width of the transition peak, and 
can consequently be determined by fitting the simulated heat 
capacity profiles to the experimentally measured one. 
The five remaining interaction parameters, , can be obtained 
indirectly from various other experiments. We have shown 
before that most anesthetics generate a freezing point depres- 
sion that is consistent with ideal solubility in the fluid phase 
and no solubility in the gel phase fl7|. Therefore we chosen 
^af = 0, and a high value for the anesthetic-gel interaction 
with LUag — ^fg ([24 1). Determination of the remaining three 
parameters [ujap, ^pg and ujpf) can, in principle, be done by 
comparing the simulated permeability (basically the average 
number of pores) to measurements. 

Even before doing this, it is possible to predict a number of 
relations, if we assume that hydrophobic matching is a major 
determinant of the nearest neighbor interactions. Firstly, if the 
added anaesthetic is strongly hydrophobic, it seems reasonable 
to assume that LOap ^ LOfg, as the cooperativity parameter, 
bjfg, is largely determined by the hydrophobic mismatch be- 
tween a fluid and gel state lipid. Similarly, this should mean 
that LOpg > Ldpf for two reasons: 1) the lipid in the gel state 
would have a larger part of the hydrophobic chains exposed to 
the water, and 2) the fluid phase is more loosely packed, mean- 
ing that defects and spaces between the lipids (i.e. pores) would 
be more likely to appear than in the tightly packed gel phase. 
Based on the considerations above, the parameter values used 
in the simulations were as follows: AH = 36400 J/mol, AS = 
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115.87 J/(mol-K), ujfg = 1326 J/mol, ujpf = 3700 J/mol, 
ujpg = 6000 J/mol, cjpa = 13260 J/mol, = 0.0 J/mol, 
and Wag = 1326 J/mol. The values for AH, AS and are 
identical to those used in f24]. As seen in the results part, 
the average number of pores form in the BLM experiments 
is in the range of 1-2 in an area of 50/im^ corresponding to 
10^ lipid molecules. In our simulations we used a matrix with 
10000 lipids that is way smaller than the experimental mem- 
brane. Therefore, for the pore formation we opted for a set of 
pore-lipid interaction parameters which on average generates a 
similar number of pores as in the experiment requiring smaller 
pore-lipid interaction parameters as one would expect in a natu- 
ral membrane. Additionally, it was necessary to disallow direct 
pore-pore interactions, so as to avoid run-away aggregation of 
pores and infinite growth of the simulation matrix. The effect 
on the results should be negligible, as typically the pore density 
in the simulation was so low as to make pore-pore interactions 
insignificant. In the experiment collision of pores is much less 
likely because of the much larger area, but may occur under 
stress when the membrane ruptures. The exact appearance of 
the pores vs. temperature profiles will, of course, depend on 
the choice of parameters. However, for a large range of interac- 
tion values, the permeability profile shows a strong peak at the 
phase transition temperature, and the position of this peak was 
found to always follow the peak of the heat capacity profile. 
Typical system sizes were 10* particles on a triangular lattice 
with periodic boundary conditions. The typical simulation was 
allowed to equilibrate for 2 • 10* to 10^ Monte Carlo cycles (lat- 
tice sweeps) two times - first without pores involved, and then 
again with pore formation allowed. The actual sampling was 
then done over another 10^ to 5 • 10^ cycles. Most of the calcu- 
lated curves presented were in addition averaged over several 
independent runs. 

3 Results 

In this paper we investigate pore formation in lipid membranes 
both experimentally with black lipid membranes (BLM) us- 
ing the Montal-Miiller technique |25J, and theoretically with 
Monte Carlo simulations. The goal is to demonstrate that pore 
formation is related to the heat capacity and the cooperative 
area fluctuations, respectively, and that anesthetics change the 
permeability in a coherent manner. 

3.1 Black lipid membranes 

In the following we show the conduction events through syn- 
thetic BLMs. We have chosen to work with lipid mixtures that 
have their transition events close to room temperature and that 
display broad heat capacity profiles. The latter makes it easier 
in experiments to adjust the temperature such that one is very 
close to the heat capacity maximum. 

Fig. [2]shows the current traces for a DOPC: DPPC=2: 1 mixture 
at 19°C (150mM KCl, pH « 6.5) for various transmembrane 
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Figure 2: Bottom: Current traces through BLMs of a 
D0PC:DPPC=2:1 mixture (ISOmM KCl, pH^6.5 at T=19°C) 
for different voltages. Top: Current-voltage relationship for the 
traces from the bottom panel. One finds a linear current-voltage 
relationship in the range from -ISOmV to +150mV. At higher 
voltages the traces become nonlinear. A total of 328 current 
traces with 30 seconds each where used for this graph. 



voltages between -220mV and H-220mV. The heat capacity pro- 
file for this lipid mixture can be found in Fig. |3]\ (bottom). 
It displays a Cp-maximum close to 19°C. In agreement with 
earlier publications (e.g., fTO'.'Sl we find quantized conduction 
events very similar to those of reported for protein ion chan- 
nels (Fig. |2] bottom) with voltage dependent amplitudes. The 
amplitudes display a linear dependence of the transmembrane 
voltage, Vm, between about -150 to + 150 mV (Fig. |2) top). At 
voltages above |Vm|= 150 mV, the conductance is found to be 
higher than expected from the linear relation. The origin of this 
non-linear current-voltage relationship is not clear but may be 
related to the influence of voltage on the phase behavior of the 
lipid membranes. It has been shown by other authors that the 
position of the melting transition can be influenced by voltage 
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Figure 3: Experiments on BLMs of a D0PC:DPPC=2:1 mixture (ISOmM KCl, pHoiS.S) in the presence (top) and the absence 
(bottom) of 15.9 mol% octanol in the membrane. A: Heat capacity profiles. The dashed line indicates experimental temperature. 
Octanol shifts the calorimetric events towards lower temperatures. In the presence of octanol the experimental temperature 
above the melting events while it is at the Cp-maximum in the absence of octanol. B: Current traces obtained at T=19°C and 
210 mV. Two representative traces are given for the two experimental conditions, respectively. Note the different amplitudes of 
the currents and the altered frequencies of current events. C: Current histograms from the curves in panel B showing on or two 
current events. D: Analysis of the frequencies of the number of pores. Symbols represent the areas of the peaks in panel C. The 
solid line represents a best fit using a Poisson-distribution. 



lISTl [32l 1331 . Thus, we suspect that the nonhnear behavior at 
higher vohages is a consequence of the influence of voltage on 
the phase transition. 

The effect of anesthetics Fig. [3]shows the influence of octanol 
on the permeation of the above membranes measured at 19°C 
and 210 mV. The top panels show heat capacity profiles (Fig. 
[3}\), representative current traces (Fig. [3jJ), current histograms 
(Fig. [3p) and pore number probability distribution (Fig. ^p) 
in the presence of 15.9 mol% octanol in the membrane (cal- 
culated using an octanol partition coefficient in the membrane 
of 200, [261). The bottom panels show the respective data in 
the absence of octanol. One can see that the presence of oc- 
tanol shifts the calorimetric profiles towards lower temperatures 
such that one is above the chain melting regime in the presence 
of octanol and in the transition regime in the absence of oc- 
tanol. The current traces display smaller current amplitudes and 
a lower frequency of current events in the presence of octanol. 
Fig. |3p shows current amplitude histograms obtained from the 
data shown in Fig. [3]3. In the absence of octanol one finds 
at least two current steps corresponding to the two maxima of 
the histogram that are different from zero. In the presence of 
15.9 mol% octanol one finds only one current step with a much 
lower amplitude. Fig. |3)D shows the areas of the current peaks 
from Fig. |3]Z! (obtained from fitting the peaks with Gaussian 
profiles). The solid lines are fits to a Poissonian distribution. 
Such distributions are expected for a low mean number of lipid 
ion channels. 

Fig. |4]shows the corresponding experiments in the presence of 



ethanol. Ethanol is a very weak anesthetic with a membrane- 
water partition coefficient of about 0.48 EtII (octanol has about 
200). Three experiments recorded at V^„=120mV are shown, 
with mol%, 9.5 mol% and 19.5 mol% ethanol in the fluid 
lipid membrane. Increasing amounts of ethanol shift the melt- 
ing profile towards lower temperatures. The corresponding cur- 
rent traces show at least two current steps at 39 and 78 pA in 
the absence of ethanol (note that the voltage is lower than in the 
octanol experiment in Fig. |3]l, while there seems to be only one 
step for 19.5 mol% of around 45pA. The current histograms 
suggests that the two currents steps seen in the absence of oc- 
tanol merge into one in the presence of ethanol. This behavior 
is different from that of octanol where the overall conductance 
decreased upon the addition of anesthetics. We did not attempt 
to make a Poisson analysis of the histograms of the ethanol 
containing membranes since the histograms in the presence of 
ethanol do not contain evenly spaced conduction levels in the 
presence of ethanol. 

In the discussion section we show that the response of the lipid 
ion channels to both octanol and ethanol is very similar to that 
reported for protein ion channels. 



3.2 Monte Carlo simulations 

In the following we attempt to rationalize the experimental find- 
ings obtained from anesthetics-containing membranes by a sta- 
tistical thermodynamics model that is describe in detail in the 
Materials section. The main idea in this simulation is to pos- 
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Figure 4: Experiments on BLMs of a D0PC:DPPC=2:1 mixture (ISOmM KCl, pH^6.5) in the presence of 19.5 mol% ethanol 
(top), 9.5 mol% (middle) and no ethanol in the fluid lipid membrane. A: Heat capacity traces. B: Current traces measured 
at 120 mV. The presence of ethanol seems to generate longer opening times. Three representative traces are given for each 
experimental condition. C: Histograms of the current events. The behavior is different as in the presence of octanol (Fig. 3] ). 
In the absence of ethanal one find two clearly visible equally spaced current steps of about 39pA that merges into one peak in 
the presence of ethanol. 



tulate a pore formation process that conserves the overall area 
of the complete membranes. The conservation of area during 
the pore formation step implies that pores can only form when 
compressing the area of the lipids. Hence, the pore formation 
process is linked to the area fluctuations of the lipid matrix. 
The area fluctuations, however, are proportional to the lateral 
compressibility. Nagle and Scott \3\ argued that the likelihood 
of finding a pore is proportional to the compressibility. Now it 
is known that the excess heat capacity of lipid melting is pro- 
portional to the compressibility |8|. We have therefore recently 
argued that the likelihood to find a pore must be intimately re- 
lated to the heat capacity [Sj. Our Monte Carlo model makes 
use of the considerations of Nagle and Scott ||3l. 

Fig. |5]shows the results of a simulated DPPC membrane. The 
simulated heat capacity profiles are compared to an experimen- 
tal profile of DPPC large unilamellar vesicles in Fig. |5]\. The 
close agreement of the two curves is not accidental but a conse- 
quence of using the parameters AH, AS and the half width of 
the experimental heat capacity profile. Fig. |5jJ shows the cal- 
culated mean number of pores in the simulation box, which is 
proportional to the conductance of the membrane. They show 
that the mean number of simulated pores reaches a maximum 
close to the transition maximum at 41°C. This behavior is ex- 
actly that found by Blicher et al. |5 1 for dye permeation through 
large unilamellar vesicle membranes and confirms the concept 
of Nagle & Scott |3| that makes use of the maximum of the 
lateral compressibility at T„i. Fig. |5p! shows representative 
traces of the pore number computed at the melting temperature. 



T„j. These traces strongly resemble those from the BLM exper- 
iments (Figs. |2]), |3]l and|4]i. Fig. |5p shows the probability dis- 
tribution of the number of pores for five different temperature 
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Figure 5: Monte Carlo simulations of a DPPC membrane. Top, 
left: Calculated heat capacity profile compared to the exper- 
imental profile of DPPC large unilamellar vesicles. Bottom, 
left: Calculated temperature dependence of the membrane per- 
meability (mean number of pores). Top, right: Four traces of 
the number of pores (corresponding to conductance) calculated 
at the melting temperature, Tm, at 41° C. Bottom, right: Prob- 
ability distribution of the number of pores. The symbols orig- 
inate from the simulation, the solid lines are fits by a Poisson- 
distribution. Simulation parameters are given in the text. 
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in this panel represent Poisson-distributions indicating that the lipid mixture. The corresponding experimental and simulated 

calculated behavior matches the experimental one (Fig. |3] bot- heat capacity profile is shown in the top panels of Fig. |6] The 

tom right). The Poissonian statistics is expected when pores maximum permeability (= mean number of pores) is found in 

form independent of each other. The lipid-pore interaction pa- the chain melting regime between the two outer Cp maxima. 
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permeability reaches a maximum in the chain melting regime. 
Simulation parameters are given in the text. 



rameters (wp/ and LOpg) are given in the Methods section and 
have been chosen such that the absolute number of pores is of 
similar order as in the BLM experiment. Since the simulated 
matrix has about 10000 lipids while the black lipid membrane 
displays a much larger area of about 1.6 • 10^" lipids, the pore- 
lipid interaction parameters would have to be about 2-3 times 
higher if the simulation matrix had the same size than the BLM 
membrane. With such parameters, however, we would hardly 
ever see a conduction event in the simulation at a matrix size of 
10000 lipids. Therefore, we have chosen smaller parameters. 
Qualitatively, however, the behavior of the pore formation rate 
with a maximum at T,„ is maintained when using the smaller 
model parameters given above. 

The above concept also works in mixtures of different lipids 
with a broad melting peak. Fig. |6] shows the mean number of 
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Figure 7; Simulations of DPPC membranes in the presence of 
anesthetics as a function of temperature and anesthetics con- 
centration. Left column: Monte Carlo snapshots in the absence 
of anesthetics. Center column: In the presence of 5 mol% anes- 
thetics. Right column: In the presence of 10 mol% anesthetics. 
The rows indicate different temperatures between 33.5° C and 
48.5 33.5° C. Gel lipids are given in red, fluid lipids in green, 
anesthetics molecules in white and pores in black. One can rec- 
ognize that the presence of anesthetics shifts the melting tem- 
perature towards lower temperatures. Simulation parameters 
are given in the text. 
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The effect of anesthetics The famous Meyer-Overton rule 1341 
states that the effectiveness of anesthetics is exactly propor- 
tional to their solubility in lipid membranes and independent 
of their chemical nature. Further, anesthetics are known to 
lower the phase transition temperature of lipid membranes. In 
a recent paper we have shown that the effect of anesthetics on 
pure lipid membranes can be well described by the well-known 

freezing-point depression law — -^jf ■ xa, where xa is the 
molar fraction of anesthetics. The derivation of this formula 
is based on only two assumption: the anesthetics molecule is 
ideally soluble in the liquid phase, and that it is insoluble in 
the solid phase. In our simulation the anesthetics molecules are 
modeled as resembling lipids that occupy lattice sites. The per- 
fect solubility of anesthetics in the fluid membrane phase can 
easily be modeled by assuming that the anesthetics-fluid inter- 
action is zero (ujaf = ), while the insolubility in the gel phase 
means that the anesthetics-gel parameter (ojag ) is very large. 

Fig. |7] shows Monte-Carlo snapshot for DPPC membranes at 
seven different temperatures evenly spaced with 2.5 degrees. 
The three columns are calculated for Omol%, 5 mol% and 10 
mol% anesthetics, respectively. Here, each anesthetics molecule 
occupies one lattice site. The melting temperature is again 
T„j = 41°C. In this figure red and green dots represent gel and 
fluid lipids, while white dots represent anesthetics molecules 
and black dots are pores. One can see how the increasing amount 
of anesthetics lowers the melting events as evident from the 
lower temperatures where one finds red gel domains. These 
simulations are further analyzed in Fig. [8] Fig. |8](top) shows 
the heat capacity profiles for the simulations in Fig. |7] They 
are found to be very similar to those found experimentally in 
the case of DPPC with octanol [17 \ and show the lowering of 
the melting temperature in the presence of anesthetics. Note, 
that the critical anesthetic dose for tadpoles corresponds to 2.6 
mol% of anesthetics in the membrane [ 17J, which is somewhat 
lower than the concentration used in this work. Fig. |8] (middle) 
shows the corresponding permeation profiles that show a sim- 
ilar shift towards lower temperatures in the presence of anes- 
thetics. Fig. [8] (bottom) shows pore formation traces for three 
different anesthetics concentrations calculated at 41 °C. These 
traces correspond to the experimental situation in Fig. |3] The 
mean number of pores is reduced. Since Monte-Carlo simu- 
lations are based on equilibrium thermodynamics the Monte- 
Carlo time scales of these traces do not correspond to a real 
time scale. This means that one cannot conclude from the sim- 
ulation that the mean life time of the pore openings in the pres- 
ence of anesthetics is longer as in the absence of anesthetics 
even though the traces seem to suggest that. 

4 Discussion 

In this paper we investigate the lipid ion channel formation in 
pure lipid membranes in the absence and the presence of anes- 
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Figure 8: Monte Carlo Simulation of DPPC membranes in the 
presence of mol%, 5 mol% and 10 mol% of anesthetics in 
the membrane. Top: Calculated heat capacity profiles show 
the shift of the melting events towards lower temperatures in 
agreement with the snapshots in Fig. 7. Middle: The calcu- 
lated number of pores ( corresponding to a permeability) follows 
the heat capacity profiles. The maximum permeability shifts to- 
wards lower temperatures upon addition of anesthetics. Bot- 
tom: Calculated traces of the number of pores indicate that 
the number of pore opening events decreases in the presence 
of anesthetics. Two traces for each anesthetics concentration 
are given. Simulation parameters are given in the text. 
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the tics. The goal was to demonstrate the origin of the pore for- 
mation from the cooperative area fluctuations. We studied pore 
formation both experimentaUy by the Montal-Miiller ||251 black 
lipid membrane technique, and by Monte-Carlo simulations. 
We show that pure lipid membranes in the absence of proteins 
display conduction events that resemble those seen in the pres- 
ence of protein ion channels. For this reason we call them lipid 
ion channels. Such events have been shown before by a number 
of aufliors |l9l[l0l[35][II][I2ISl[l3l|37l[Ml[39l[Ilia. Alfliough 
there are obviously many reports about such events it is surpris- 
ing that such behavior is not better known. While the pores in 
lipid membranes are shown to exist in the absence of proteins, 
the conduction of ions through proteins cannot be measured in 
the absence of lipid membranes. Since the events are so similar 
to the protein ion channel recordings one requires an unam- 
biguous criterion for how to distinguish the events originating 
from proteins and lipid membranes. To our knowledge, no such 
criterion exists. Possibly, site-specific mutations of proteins re- 
lated to changes in channel activity can be considered as an 
indirect criterion, but this may not be satisfactory because this 
only demonstrates that the protein are somewhat involved in the 
conduction process (cf. last paragraph of discussion). 
The importance of phase transitions in lipid membranes for the 
permeability has been pointed out before Il2l|3][l0ll2l|4ll21. The 
cooperative fluctuations in area render the pore formation more 
likely because the lateral compressibility is higher. This effect 
resembles the phenomenon of critical opalescence described by 
Einstein f40l. He calculated the work to move one molecule 
from one part of the volume to another in in order to obtain the 
density fluctuations. In a similar manner, the increased fluctu- 
ations in the 2D membrane render it easier to move lipids such 
that pores can be created. In the present paper we used this con- 
cept to describe the pore formation process by a simple statisti- 
cal thermodynamics model in which lipids are moved while si- 
multaneously reducing the area of the surrounding lipids. Since 
the change of the intensive thermodynamic variables (e.g., pres- 
sure, lateral pressure, chemical potentials of protons, of cal- 
cium or of anesthetic drugs) influences the melting behavior 
of membrane and thus the fluctuations in the membrane, it is 
straight forward to explore how changes in these variables in- 
fluence permeability. For instance, it has been shown that volt- 
age influences the phase transition temperature ||32l |33l . As 
evident in our measurements, the mean conductance through 
pores in lipid mixtures is linear in the range between about - 
150 to H-150mV. For \Vm\ > l50mV the current-voltage rela- 
tion becomes non-linear indicating that high voltage influences 
the state of the membrane and induces more conduction events. 
It has further been shown that changes in calcium concentra- 
tion f36] [T3]| and pH changes fTTl influence to occurrence 
of lipid ion channels. Here, we have used the chemical poten- 
tial of anesthetics in the membrane as a variable. Our present 
work is continuation of a recent publication from our group fS) 
that shows how anesthetics influence the permeation process by 
influencing the melting process. Here, we have demonstrated 



this for the case of octanol and ethanol, which both shift tran- 
sitions to a similar degree when the membrane concentration is 
the same. Both drugs lead to a reduction of the overall conduc- 
tance. 

Our simulations described the main feature found in experi- 
ments: The relation of the permeability and the heat capacity 
leading to maximum permeability at the heat capacity maxi- 
mum. This confirms the main assumption of our study, which 
is that the cooperative area fluctuations are responsible for the 
occurrence of the lipid ion channels. It further describes quan- 
titatively the effect of anesthetics on the heat capacity of the 
membrane. The anesthetics following the Meyer-Overton cor- 
relation also generate freezing-point depression. The simula- 
tion accurately describes this and therefore also describes the 
influence of anesthetics on the permeation. However, our pore 
formation model made the simplifying assumption that a pore 
has the size of one lipid. This is a reasonable approximation 
since we found in a recent paper that the conductance through 
the membrane is well described by assuming an aqueous pore 
of a diameter of 0.7nm, which is exactly the size of one lipid. 
The measurements in the presence of anesthetics show, how- 
ever, that the pore size itself can vary and that it is different in 
the presence and absence of octanol and ethanol. This effect is 
not contained in our simulations. 

The present Monte Carlo simulation is a coarse grain statis- 
tical thermodynamics model. Such models have successfully 
been used to describe cooperative phenomena in membranes 
HSl [m EOl ED ED EH. Typically, the number of states for 
each lipid is reduced to two (Ising model) or 10 (Pink model). 
Naturally, this is only a very approximate description on the 
molecular scale. The reason why the lattice models work so 
nicely is that they are designed to describe cooperative behav- 
ior of many molecules - in our case the formation of domains 
and the cooperative fluctuations. Microscopic details of single 
lipids and pores described here by lattice sites are not expected 
to be accurate - and also to be of small relevance. As men- 
tioned above, our simulation is designed such that a pore has 
always the size of one lattice site. The experiment on octanol 
shows that the pores in the presence of anesthetics are smaller 
than in the absence (currents at identical voltage are « 90pA 
in the absence of octanol and only about «15pA in the pres- 
ence of 15.9 mol% octanol). Qualitatively, however, the sim- 
ulation describes the likelihood of pore formation and thereby 
permeation quite well. In a recent paper we have shown that 
the permeation of fluorescence dyes through lipid vesicles is 
within experimental accuracy proportional to the heat capacity 
|5|. This paper also shows the quantized conduction events in 
BLMs that could be abolished by addition of octanol. Our sim- 
ulations reproduce this close relation between heat capacity and 
conductance - and is consistent with the quantized nature of the 
conduction process. 

It is likely that many other variable changes would influence 
the occurrence of lipid ion channels, e.g., proteins that asso- 
ciate or insert into membranes. It is known that band 3 pro- 
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tein of erythrocytes increases the mehing temperature of model 
membranes Ii4rj . Cytochrome bs, on the other hand lowers 
transitions ll42ll . Cytochrome c binds peripherally and increases 
the transition temperature (21] . Depending onwhether the ex- 
perimental temperature is found below or above the melting 
temperature, one would expect that the respective proteins can 
either increase or lower the likelihood of lipid ion channel for- 
mation even if the proteins are not known to form ion channels. 
The mechanism of anesthesia has been under much debate and 
is unexplained until today. Until the mid 1970's lipid models 
were most popular because they can be easily related to the 
famous Meyer-Overton rule ll34ll that states that the effective- 
ness of anesthetics and their partition coefficient are linearly 
related, independent of the chemical nature of the anesthetics 
molecules. Ever since the development of the patch clamp 
methods and the study of protein ion channel events, a lot of 
attention has been on the binding of anesthetics to receptors. 
This is due to the finding that ion channel events that have been 
attributed to proteins are influenced by anesthetics. Unfortu- 
nately, such interaction with proteins is not generally in agree- 
ment with the Meyer-Overton correlation. Further, we have 
shown here and in fSl that the finding of quantized events does 
not prove the action of a protein. Further, the influence of anes- 
thetics on the conduction events does not prove a specific bind- 
ing of the anesthetic drug to a receptor We have shown here 
that octanol and ethanol influence the channel events even in 
the complete absence of proteins. 

In this context it is tempting to compare the effect of alkanols 
on protein ion channels with our findings. The influence of oc- 
tanol on current fluctuations of lipid ion channels can be com- 
pared with those of neuronal nicotinic acetylcholine receptor 
(AChR, see, e.g., [43]), which has been used for many years 
as a model to study the influence of alcohols ll44l l45l l46l l47l 
l48l l49l l50l . It has been shown that neuronal activity of nico- 
tinic acetylcholine receptor (AChR) induced channels under 
the influence of octanol is markedly inhibited ||43| . which was 
observed as a decrease in current fluctuation amplitude from 
those channels. Ethanol, on the other hand was found not to be 
neutral for many protein ligand-gated ion channels such as N- 
methyl-D-aspartate (NMDA), serotonin (5-HT(3)), glycine and 
GABA receptors |51 52|. The inhibitory action of ethanol was 
also observed in voltage-gated Ca2H- channels li53l and very re- 
cently also in potassium channels, which were seen as very im- 
portant targets of ethanol |54|. Similar findings were reported 
for experiments with AChR channels subjected to lOOmM of 
ethanol 1551 . The above reports demonstrate that the effect 
of octanol and ethanol on membrane protein receptors is quite 
similar to that found by us for the effect of these alcohols on 
lipid ion channels. In this context it is also interesting to note 
that protein ion channels have been shown to respond to the 
phase behavior of lipids Il56l or to different lipid environments 

El. 
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